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Abstract. The transport properties of quantum dots with up to N = 7 electrons 
ranging from the weak to the strong interacting regime are investigated via the 
projected Hartrce-Fock technique. As interactions increase radial order develops 
in the dot, with the formation of ring and centered-ring structures. Subsequently, 
angular correlations appear, signalling the formation of a Wigner molecule state. 
We show striking signatures of the emergence of Wigner molecules, detected in 
transport. In the linear regime, conductance is exponentially suppressed as the 
interaction strength grows. A further suppression is observed when centered-ring 
structures develop, or peculiar spin textures appear. In the nonlinear regime, the 
formation of molecular states may even lead to a conductance enhancement. 
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1. Introduction 

Semiconductor quantum dots (QDs), frequently referred to as artificial atoms, are 
nanometer-sized structures whose conduction electrons are confined in all the three 
spatial dimensions [TJ[21[3]. In these systems a two-dimensional electron gas, formed 
at the interface of a heterojunction, is depleted by chemical etching or electrostatic 
potentials in order to form an isolated region, connected to external reservoirs by tun- 
nel barriers. For a small number of particles N, the potential can often be considered 
as harmonic [TJ [U [3] . 

In analogy to atomic systems, quantum dots can be probed optically by studying their 
absorption or emission spectrum [JJ . Additionally, the study of transport properties 
is a source of information for quantum dots embedded into an electronic circuit [TJ. 
The current flow proceeds by tunnelling events once a bias voltage is applied to the 
external reservoirs and the presence of an external gate voltage allows to tune the 
number of excess electrons in the dot with respect to a neutral configuration. 

Theoretically, the study of correlated quantum dot states is a challenging many-body 
problem: even the fairly simple case of N = 2 can be solved exactly only in specific 
regimes [5], whereas in general one has to resort to semi-analytic methods [6] or 
approximate WFs [7j. For N > 2, several numerical methods have been employed. In 
increasing order of computational complexity they range from unrestricted Hartree- 
Fock methods | | [1 E] and density functional theory PJ [TJ HI UJ3 EE], 
to projected Hartree-Fock (PHF) [171 [TBI El EH HOI HH [22 [23], random phase 
approximation [24j [25], quantum Monte Carlo [TBI EH EH EH [29], and exact 
diagonalization [SBElliaiMlElliaiiaEa^ Recently, PHF techniques 

have been used by our group for the study of correlated quantum dots. While 
retaining the flexibility of an unrestricted Hartree-Fock approach, PHF allows to 
overcome the limitations due to symmetry broken solutions and to efficiently obtain 
dot wavefunctions (WFs) with the correct spin and angular momentum, showing 
correlations beyond the mean field level j21j [22] . 

Circular two-dimensional systems realized e.g. by pillar quantum dots display 
interesting features depending on the ratio A between the strength of the typical 
interaction and the strength of the confining potential - see Sec. 12.11 for a precise 
definition of A. When A is small, atomic-like effects due to quantum mechanical 
confinement such as the formation of shell structures have been observed [40] [UJ [42] 
and explained at the mean field level p] 143] - For increasing A, several numerical 
investigations have shown the emergence of correlated electron states and the 
occurrence of Wigner molecular states gllSllMlMllSTllSHllSSliaililSllSlIIl 
[4l[TJ[43[4![4S[5rj[33[^ the finite-size 

analogue of Wigner crystals [59j [60] , characterized by correlations beyond the mean 
field. 

The transition towards a molecular state occurs smoothly [26J [TS] [27l [28] [29] [47l [58] . As 
interactions increase, the dot WFs cross over from weakly correlated states at small 
A to Wigner molecular states, characterized by strong correlations. The transition 
occurs in two phases: 

• At small A, correlations begin to develop and ring-like structures develop in the 
dot WF. For TV > 6, also centered structures with one or more electrons in the 
dot center may form. 
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• For higher A, angular correlations begin to appear as the dot enters the incipient 
Wigner molecule regime |27[ I28| . Increasing A further, the dot WF represents a 
rotating Wigner molecule |45j and the electrons localize around the equilibrium 
positions of a classical Coulomb molecule [HI |4"9"] . 
Correlations are particularly relevant also in one-dimensional systems which display 
an analogous transition towards the Wigner molecule [5TJ [HH [531 [Ml ES] ■ 
The experimental observation of strongly correlated states in quantum dots has 
attracted considerable interest. In pillar quantum dots, inelastic light scattering 
experiments have shown signatures of correlated quantum states [571 166) . Scanning 
tunnelling spectroscopy experiments for the imaging of correlated quantum dot WFs 
were recently performed and theoretically analyzed [551 1561 I67| . Also transport 
properties can yield information about correlated states. In one dimension, the 
influence of correlations on the transport properties is predicted to be particularly 
important [63 [69l EOj. Recently, experimental evidence of the formation of few- 
electrons Wigner molecules has been reported in carbon nanotubes [71] . 
Also spin correlations can heavily influence the transport properties of quantum dots, 
even in the absence of an applied magnetic field. In quantum dots with asym- 
metric tunnel barriers, the degeneracy of spin multiplcts may lead to asymmetric 
current-voltage characteristics [72] . Another notable example is the type-II spin block- 
ade [731 174] , which occurs in the linear transport regime when the absolute value of 
the difference between the total spin of initial and final dot states exceeds 1/2 and 
leads to zero sequential current through the dot. 

In this paper we investigate the transport properties of quantum dots in the presence 
of strong correlations. In such a regime, a mean field treatment in the spirit of the 
so called "constant interaction model" [40l [75] is clearly not viable. Indeed, one has 
to resort to more precise techniques to obtain the spectrum and the WFs. Numerical 
studies of the transport properties, similar to the one proposed here have been per- 
formed in the past employing exact diagonalizations for N < 3 electrons in circular 
QDs [76l [771 EH] and N < 4 electrons in one-dimensional quantum dots [79] . These 
works, however, were not focused on the signatures due to Wigner molecules in the 
transport properties. 

In the present work we numerically investigate the transport properties of pillar quan- 
tum dots beyond the constant interaction model. Our model is that of N interacting 
electrons confined to a two-dimensional plane and further subject to an in-plane har- 
monic potential. More refined models, including effects due to a finite thickness of 
the dot and to heavy doping in the reservoirs, have been recently proposed |80j . In 
this work we will neglect such effects, addressing systems in which the screening is 
moderate (a strong screening may hinder the formation of Wigner molecules |61[l81j). 
We use the PHF method in order to estimate the correlated dot WFs for 4 < TV < 7 
in a range of A which allows to observe the transition between liquid-like and molecu- 
lar electron states. Sequential tunnelling rates are numerically evaluated and the dot 
conductance is obtained using a rate equation. 

Our task is to understand whether or not peculiar signatures in the transport proper- 
ties may be detected as a consequence of the transition towards the Wigner molecule. 
According to the results presented in this paper, the answer is affirmative. 
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In the linear transport regime, qualitative modifications of the dot ground state (GS) 
WFs induce a peculiar suppression of the conductance. Such qualitative modifica- 
tions may be induced either by the formation of centered ring-like structures or by the 
emergence of peculiar spin patterns in the dot WF. Both cases are presented in this 
paper. 

Signatures of the transition can also be seen in the nonlinear transport regime. We 
have found that the tunnelling rate through an excited state of the dot may be in- 
creased strongly by the formation of a Wigner molecule. 

The features described above are genuine hallmarks of the formation of Wigner 
molecules in the dot and can be expected to be observable in experiments. 

The outline of the paper is as follows. In Sec. [2] we introduce the model and the PHF 
method, we discuss the tunnelling Hamiltonian and the rate equation for calculating 
the current. Results are presented in Sec. EH Here, after discussing in detail the 
occurrence of Wigner molecules, we show results for the conductance in both the linear 
and nonlinear regimes. Conclusions are presented in Sec. [4] |Appcndix A| contains the 
derivation of the tunnelling Hamiltonian while the dot tunnelling rates are discussed 
in |Appendix B| 

2. Model and methods 

2.1. Quantum dot 

In a pillar quantum dot [40j EH 021 [43] electrons are confined to a thin disk of 
semiconducting material, represented by the red region in figure [TJ The dot is 
embedded between tunnel barriers located around z = z~e,zc, with z the axial 
direction. The tunnel barriers couple the dot to the external emitter and collector 
contacts [2]- A metallic gate is assumed to surround the dot region (not shown in 
figure [1]) and allows to shift the dot energy levels as a suitable gate voltage V g is 
applied to it. Due to the strong confinement along z, the motion of electrons is 
essentially restricted to the (x, y) plane. Electrons are further subject to a lateral 
confining potential with rotational symmetry around the z axis, as appropriate for 
the cylindrical quantum dots studied in this paper. For small dots, containing few 
electrons, this potential is well approximated by a parabolic one [2]. The Hamiltonian 
for A" interacting electrons is = with (h = 1, boldface denotes vectors) 

AT r , * 9 i 2 N-l N 1 

h l 2m * 2 J ^Sr {-{ IB, - R,| 

where = (xi,yi) is the i-th electron coordinate and Pi its momentum. Here, 
— e and m* are the electron charge and effective mass respectively. Furthermore, u) 
is the confinement energy, e (£r) the vacuum (relative) dielectric constant. We will 
consider a bare Coulomb potential for the interaction term H^' , neglecting both finite- 
thickness effects and screening due to heavily doped contacts. Such effects modify the 
interaction potential producing deviations from the behaviour for both short and 
long inter-electron distances r [501. Finite thickness effects would also produce a 
renormalization of the gate voltage V g [50] . Our calculations are therefore valid for 
systems characterized by weak screening. 
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Figure 1. Schematic representation of a pillar quantum dot. The red thin disc 
is the quantum dot (D), connected via tunneling barriers at z as z-q,zq (yellow 
parts) to emitter (E) and collector (C) leads, represented in blue. 



Expressing lengths in units of £q = (to*w) -1 / 2 and energies in units u, the Hamiltonian 
becomes 

N r~9 -91 Tr N-l N 



(0) 



E 

i=l 



p r 



w 



,(2) 



with f j = Hi/£o and p.j = P 
A 



i*0 



The dimensionless parameter 



= — (3) 

measures the Coulomb interaction strength. It is the ratio between the effective 
length scale £o and the effective Bohr radius a,g = 4irsoS Y /m*e 2 . Experimentally, 
the interaction strength A can be modified by tuning the confinement strength to via 
electrostatic gates. 

In the rest of the paper, we will concentrate on GaAs quantum dots, where e r = 12.4 
and m* = 0.067m e with m c = 9.1 ■ 10 -31 kg. In this case, expressing u in meV, one 
has A w 3.46 VmeV /y/u>. Weak (strong) interactions occur for A < 1 (A > 1). In the 
absence of interactions (A = 0) the problem can be solved exactly [52] ■ The eigenstates 
of Hj^ are Fock-Darwin (FD) states labelled by a principal quantum number n > 0, 
by the electron angular momentum (z component) I G Z and by the electron spin z 
component s z = ±1/2. The corresponding WFs are denoted by f n ,i,s z (r) and the spin 
degenerate energy spectrum is given by E ny i^ s _ = w(2n + \l\ + 1). In the presence of 
interactions, the problem cannot be tackled analytically if N > 2 and one has to use 
numerical techniques. It is important to notice the symmetries of H^: it commutes 
with the total angular momentum (z component) L, the total spin S and the total 
spin z component S z . As a consequence, their eigenvalues can be used to label the dot 
energy spectrum and WFs. These are obtained by means of the PHF technique which 
has been extensively described in [3TJ[22]. Here, we briefly outline the procedure. For 
a given particle number N and each value of —N/2 < S z < N/2, the dot WFs are 
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first approximated as single Slater determinants \N, S z ) made up of iVj (-/VjJ orbital 
with spin s z = 1/2 (s z = -1/2) where S z = (N^ - N{)/2 and N = iV T + Ny 
Orbitals are variationally optimized with the spin and spatially unrestricted Hartree- 
Fock method |5] which produces several stationary states \N,S z )i, in general neither 
cigenstates of L, nor of S 2 . Projection operators Pl,s are subsequently applied to 
\N,S z )i to restore the symmetries broken due to the single Slater determinant ansatz. 
As a result, correlated WFs 

\N,L,S,S z ) i =P LtS \N,S z ) i (4) 

are obtained. The state in ^ cannot be represented as a single Slater determinant 
and contains correlations beyond mean held. The dot ground state is obtained as the 
state which minimizes the energy 



P . J jW S z \HP LtS \N,S z )j \ 

M *° ,..<< . ,\ .(^S^Pl^S^ J 

Here, Lq, So, S z q is t ne set of quantum numbers which minimize ([5]) and which label 
the dot ground state for N electrons, with WF \N, Lq, So, S z q). In a similar fashion 
one defines excited states within the PHF method. For instance, the first excited state 
is given by 

E NM , = ^ i{NiStlhMNtSa)i > Z» Mo <j , (6) 

with L\, S%, and S z \ determined by the minimization procedure. As a consequence 
of the correlations introduced by the projection technique, energies lower than those 
obtained by unrestricted Hartree-Fock are achieved [HI [55] . 

2.2. Tunnelling rates 

As shown in |Appcndix A[ assuming cylindrical symmetry about the z axis and 
separability of longitudinal and transverse motions, the tunnelling Hamiltonian 
between the quantum dot and the leads is [83] 

A *= E E^ 6 U^+h.c (7) 

where a = E (a = C) for the emitter (collector), £ a and r\ collectively denote a 
set of single particle quantum numbers for the lead a and the dot respectively and 
Tg°^ is the tunnelling amplitude. In the following, we will choose the FD states 

r) = |?1t; D \ li^ , Sz))' | as a basis of single particle states for the dot. In the case of a 

pillar quantum dot one obtains r^ aiV w t^ a '8 VatV , see (|A.7|) . The choice of FD states 
is not restrictive: indeed every orthonormal and complete basis for the single particle 
states of the dot produces identical results, as shown in |Appcndix C| The fermionic 
operator for the lead a is c a ^ a , while d v is the one for the dot. Leads are treated as 
nonintcracting Fermi gases with the Hamiltonian 

H a = J2 E ^ a )^U a C a , ia , (8) 
with energy spectrum E a (£ a ). 

Our task is to evaluate the sequential tunnelling rates between initial (|2b)) and final 
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(I^T)}) dot states with energies Ej D and Ef , respectively. As shown in | Appendix B| 
the rates are obtained by tracing out the degrees of freedom of the leads, and have 
the general form 

rzo^o = E E rgj^ , 0) 

o=E,Cp=il 

where p = +1 (p = —1) represents tunnelling into (out from) the dot via lead a. They 
are 

4^ D =r (a) |Op| 2 / P (MD-/ia), (io) 

where = 27rD Q |^ a - ) | 2 is the bare tunnelling rate with V a the density of states 
of lead a and f p (E) = pf{E) + (1 - p)/2, where f(E) = [1 + exp^E)]" 1 is the 
Fermi distribution at inverse temperature /3 = 1/fceT (fce the Boltzmann constant). 
The chemical potential of the dot is fi\-y = Ejr u — Ex D and those for the leads are 
Ma = Mo + Sfi a . Here, 5[A a is a shift due to the presence of a bias voltage V . In the 
following, symmetric voltage drops will be assumed at the barriers, with SfiE = eV/2 
and 5[ic — —eV/2. 

Interaction effects are embodied into the term 

Ox = E^dI^IXd) ; O-i = E(^d|4I^d) , (11) 

which can be evaluated numerically once the initial and the final dot states have been 
obtained by means of PHF. For A = 0, one can only have \O p \ 2 = 0, 1 depending on 
the initial and final dot states. For A > 0, on the other hand, \O p \ 2 is not limited to 
these two extreme cases. Note that \O p \ 2 contains interference effects between differ- 
ent FD orbitals. 



2.3. Rate equation 

Using the tunnelling rates one can set up a rate equation for the occupation 
probabilities P% of the dot states \T) (in this section, we omit the subscript D for 
simplicity) 

^=Y,M xjPj . (12) 
J 

The rate equation is a powerful and standard tool to study the transport properties 
of quantum dots, especially in the sequential regime [75] , 
The transition matrix Mxj is defined as 

M XJ = Yj^x (1 - S Nx , N:r ) + Wj^x5n^ Nj Hl^J 
Mxx = - E Mx ' x ' 

X'jtX 

with r j—,x given in ([9]) . In order to take into account dissipation effects on the excited 
states, we introduced a phenomenological relaxation rate 

W x ->j = W if Ej < E x ; 

Wx^j = We~ p{Bj - Bl) if Ej>E T . (13) 

In the stationary regime, the left hand side of (JT2J) vanishes and the rate equation 
reduces to a standard linear system of equations for the stationary occupation 
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probabilities of dot states P% which can be easily solved by means of singular value 
decomposition since dct(M-jj) = 0. The solution is uniquely determined by imposing 
the normalization condition Y] T Px = 1. Once the dot occupation probabilities are 
obtained, the stationary current I^ a ' through barrier a can be calculated with the aid 
of the barrier-resolved tunnelling rates (p~0|) as 

^'-EEE^. d4) 

In the stationary regime, /( E ) = — /( c ) = /. The differential conductance is defined 
as Q = dl/dV. 

2.4- Quasiparticle wavefunction 

Useful information about the dot states can also be extracted from the quasiparticle 
WF (QPWF) [56] 

<p{T)= (*bl#»|Zb) (15) 

s z =±l/2 

where 

n>0,l 

is the dot field operator, with f n ,i,s z { r ) t ne FD WFs, and the final dot state \Ft>) has 
one extra electron with respect to the initial one |2"d)- The squared modulus |v(r)| 2 
is proportional to the probability density of tunnelling into the dot at position r. 
The QPWF is the analog of the single particle WF of a tunnelling electron for the 
case of an interacting quantum dot: for A = it simply reduces to the WF of the FD 
orbital occupied by the tunnelling electron. 

For a transition from the state |7V, L, S, S x ) to \N + 1, L', S', S' z ), the QPWF has the 
general form in polar coordinates r — > (r, 9) [56] 

v(r,9)=e WAL \cp(r)\, (17) 

where AL = L' — L. Furthermore, |v?(r)| cx r' Ai ' if r — > and |y(j")| — > for r — > oo. 

3. Results 

In this section we present results for a GaAs-based quantum dot with 4 < N < 7, 
with parameters e r = 12.4 and m* = 0.067m e where m c = 9.1 • 10 -31 kg. 
For the PHF calculations, we use a truncated basis consisting of 75 FD states per 
spin direction. Projection operators are numerically implemented with a fast Fourier 
transform over 256 samples. For further details, see [22]. The ground state and first 
few excited states are obtained for interaction strengths in the range 1 < A < 2.8. 

3.1. Molecular states of electrons 

In table [T] the quantum numbers of the many-body ground states of the dot for 
increasing values of A are shown as derived from the PHF procedure. Dot states 
consist of multiplets, degenerate on S z and on L = ±Lq. Dot quantum numbers are 
constant throughout the whole range of interaction strengths 1 < A < 2.8 considered 
in this paper. They agree with the results of more refined exact diagonalizations [39] . 
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Table 1. Quantum numbers of the dot GSs as a function of TV in the range 
1 < A < 2.8. 



N 


L 


S 


s z 


4 





1 


0,±1 


5 


+ 1 


1/2 


±1/2 


6 











7 


±2 


1/2 


±1/2 



As the interaction strength increases, the dot GS WFs undergo profound 
modifications, crossing over from weakly correlated states at small A to Wigner 
molecular states, characterized by strong correlations among electrons, for higher 
A [261 H2 EH Ell EH [47l El]. The crossover is smooth and occurs through two 
phases [27l EH [29] . 

First, radial correlations begin to develop. As a result, ring-like structures are formed. 
In addition, for TV > 6 also centered structures appear, with the localization of one or 
more electrons in the center. 

As A is increased, angular correlations begin to develop entering the incipient Wigner 
molecule regime [57J EH] ■ Eventually, for strong interactions, the dot WF becomes a 
rotating Wigner molecule, with electrons localized around positions corresponding to 
those of classical charged particles parabolically confined [H] 25] . Such states are the 
analogue of the Wigner crystal [55J but with a finite size. Angular correlations cannot 
be resolved in a rotationally invariant system but can be characterized by studying two- 
body angular correlation functions, which show an oscillatory behaviour [27 ] [28 ] [29], 
Also the WFs calculated with PHF show a behaviour in qualitative agreement to the 
above results. Let us begin to introduce the spin-resolved one-body electron density 
p\- (r), defined for a normalized dot state | TV, L, S, S z ) as 

ft (r) =(N,L, S, S z \i>l (r)&, (r) |TV, L,S,S Z ) (18) 

and ij) Sz (r) defined in (fi"6]) . In this section we want to illustrate the most relevant 
aspects of the transition towards the Wigner molecule and will not enter into details 
about the spin structure of such states. Therefore, we only need to consider the total 
charge density, summed over the spin: 

p1 (t)= J2 p'l'W- ( 19 ) 

s z =±l/2 

Since the dot WFs are eigenstates of the angular momentum, introducing polar 
coordinates r — > (r, 9) one has p\ (r) = p x (r) . A plot of p\ (r) for different values 
of A is represented in figure [5] For TV = 4 and TV = 5, with increasing interactions 
the density is depleted in the core of the dot and a sharp ridge is formed at its edge, 
suggesting the formation of a ring-like structure. The position of such ridge moves 
outwards as the interaction strength increases. Also for TV = 6, 7 a ridge develops at 
the edge and moves outwards for increasing A. Additionally, for TV = 6 (TV = 7) the 
density develops a bump for r k when A > 2 (A > 1.8). This behaviour is consistent 
with the formation of a centered ring structure. As we shall see in the next section, 
this latter rearrangement of the WF produces detectable signatures in the transport 
properties. All these findings show that, for increasing A, radial correlations among 
electrons get more pronounced. 




Figure 2. One-body density pi(r) (units: £g 2 ) of the dot GS as a function of 
the distance from the dot center r, for different interaction strengths A (shown 
in the inset) and number of particles N: (a) N = 4; (b) N = 5; (c) N = 6; (d) 
N = 7. For dot parameters, see text. As discussed in Sec. 12.11 here and in the 
following r is normalized to £o- 



In order to investigate the development of angular correlations and the emergence of 
a Wigncr molecular state, one can introduce the two-body correlation function 

pi ' = (N,L,S,SSt(r)4>ly)^M)^ Sz (r)\N,L,S,S z ), (20) 

connected to the pair distribution function g Sz ,s' (r, r') by p s 2 z,Sz (y]y') = 

p{* (r)pl* (r')g SziS > z (r, r') [60]. It is proportional to the conditional probability of finding 
one electron with spin s z at r. provided that another electron with spin s' z is at r'. For 
the qualitative discussion in this section, we consider the total two-body correlation 
function 

P2 (r,r')= J2 E P^'M- (21) 

«*=±l/2 s' z =±l/2 

An example of the spin structure of the Wigncr molecules will be discussed by 
employing p2 Z,s '(v) in Sec. I3.2.T1 in connection with transport results. 
A natural choice for studying p 2 (r,r') is to fix r' at one point on the ridge of the 
one-body density: |r'| = r (N 7 X) and 8' — 0, where r (N 7 X) denotes the position of 
the off-center maximum of pi(r) for N electrons at interaction strength A. 

Figure [3] shows p2{&) = P2(ro(N, A), 8, ro(N, A), 0) as a function of 8 for different 
values of A. For weak interactions (A ~ 1, red and green curves), the correlation 
function is almost flat except for the "Fermi hole" at 8 = 0, 2ir, essentially induced 
by the Pauli exclusion principle. This confirms that correlations among the electrons 
within the ring are weak. For increasing A, the depletion at 8 = 0, 2ir gets more 
pronounced, signalling the increased importance of dynamical correlations. Even 
more important, at the highest values of A considered, P2(&) develops an oscillating 
structure, consisting of JV— 1 maxima for states with TV = 4, 5 and with N — 2 maxima 
for N = 6, 7 electrons. This is is consistent with the discussion above, namely that 
angular correlations "lag behind" and appear for values of A higher than those at 
which radial correlations get sizeable. Combining the information gathered from the 
electron density and the two-body angular correlation function, one can expect that 
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2?r B 2-k 

Figure 3. Angular dependence of the two-body density functions p2(S) = 
p2(ro(N, X),e,r (N, A),0) (units: £g 4 ) for the dot GSs for different values of 
A and (a) N = 4; (b) N = 5; (c) N = 6; (d) N = 7. Color-codes for A are the 
same as in figure [2] 



for strong interactions the dot WFs for TV = 4 and TV = 5 have the structure of a 
square and a pentagon, respectively. For N ~ 6 and TV = 7, they resemble a centered 
pentagon and a centered hexagon. 

This is confirmed by figure HI which shows a density plot of P2(r, r') for the dot 
GSs as a function of r in the (x,y) plane. The white cross denotes the position of 
r', which is the same as in figure El Around r', the presence of the Fermi hole is 
clear. Strong radial and angular correlations are observed, confirming the structures 
for TV = 4 (square), TV = 5 (pentagon), TV = 6 (centered pentagon) and TV = 7 
(centered hexagon). For comparison, density plots for A = 1.2 are shown in figure [5] 
For such a smaller interaction, radial and angular correlations beyond the Fermi hole 
are undetectable. The situation is more reminiscent of a liquid-like behaviour. 
Even though the discussion has been focused on the dot ground state, also the WFs 
for the excited states behave in a similar manner. 

The above results confirm that PHF is able to capture at the qualitative level all 
the relevant correlations of the dot WFs and to produce Wigner molecular states. In 
this respect, we note that the onset in A for the development of strong radial and 
angular correlations in the WFs predicted by the PHF method seems to be smaller 
than the one found with other techniques. As an example, for TV = 6, both exact 
diagonalization [56] and density functional [12] calculations predict the localization of 
one electron in the dot center for A « 8 while from the PHF calculations one would 
obtain A ks 2. A similar tendency to underestimate the crossover in A for the transition 
between different dot GSs has already been observed in earlier studies of PHF [22\. 
Since the qualitative changes of the dot WF are correctly captured by PHF, we expect 
that the transport results described below will be at least qualitatively correct. 
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Figure 4. Density plot of the two-body density p2(r,r') (arbitrary units, 
topographic color-scale: red high and blue low values) for the dot GS with (a) 
N = 4; (b) N = 5; (c) N = 6; (d) N = 7 and A = 2.8. The coordinate r' of the 
probe electron is chosen on the external ridge of pi (r) and is marked with a white 
cross in each panel. See text for other physical parameters. 



3. 2. Transport properties 

In this section we will show how modifications of the WF, occurring in the transition 
from a liquid to a molecular character, can be detected in the transport properties. In 
the rest of the paper, we assume symmetric tunnelling barriers, with r( E ) = r' c ' = Tq 
in (flC)]) . Typical values for Tq are of the order of some MHz. 



3.2.1. Linear transport We start considering the linear regime (V —> 0), which 
provides information on the dot ground states. A plot of the linear conductance 
Q as a function of V g , calculated solving numerically (fl"2"|) in the stationary regime is 
shown in figure [Ha). It has been calculated for A = 1.2. The conductance exhibits 
the well known Coulomb oscillations: conductance peaks are separated by regions 
where the dot is in the Coulomb blockade regime and transport is forbidden Q3 [75] . 
Peaks occur when the chemical potential of the dot is aligned with the electrochemical 
potential of the leads fj,r> = /Act) which is satisfied for a given transition TV <-> N + 1 
by suitably tuning V s . Since no simply induces a constant shift of the position of the 
linear conductance peaks in V s , we assume = 0- Turning to stronger interactions 
A = 2.4, figure [6jb), the linear conductance decreases. The observed suppression of 
Q as A is increased can be interpreted as due to the increased difficulty to tunnel 
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Figure 5. Same as in figure|4]but with A = 1.2 in all panels. 



into (or out from) an electronic system with strong Coulomb repulsion. However, the 
conductance peaks for the transition 5 <-> 6 and 6^7 have been suppressed much 
more than that corresponding to 4 <-> 5. 

In order to investigate this behaviour more systematically, the heights of the conduc- 
tance peaks are shown in logarithmic scale as a function of A in figure E^c,d). For 
the transition 4 <-> 5 (circles), a single slope is observed, signalling an exponential 
suppression of the conductance as A increases. On the other hand, for 5 <-> 6 (squares) 
a bimodal behaviour occurs, with a slope for A < 2 that is very similar to the one 
found for 4 <-> 5. A steeper slope is found for A > 2. The conductance peak for the 
transition 6 <-> 7, sec figure EJd), shows a behaviour similar to 5 <-> 6: a smaller slope 
for A < 1.6 and a steeper one for A > 2. 

In order to interpret these behaviours we can deduce more precise information about 
the tunnelling of electrons from the QPWF, sec (fl~5)) . Figure [7] shows its modulus 
\<p(r)\ for the transition between dot GSs N —> N + 1 with N = 4 (a), N = 5 (b) 
and N = 6 (c) and increasing values of A. Since all these transitions have |AL| 7^ 
(sec table Q}, the WF exhibits an off-center maximum and is small around the origin, 
hence tunnelling is strongly suppressed in the center while it is enhanced at the edge 
of the dot. The above transport results are now explained by considering both the 
shape of the QPWF and the structure of the WF of the dot GS for each N, discussed 
in Sec. O 

On the one hand, by comparing the WFs for two subsequent dot GSs one can estimate 
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Figure 6. (a) Plot of the linear conductance Q as a function of V g (units mV) 
for T = 500 mK and A = 1.2. (b) Same as in (a) but for A = 2.4. (c) Plot of the 
conductance maxima (log scale) as a function of the interaction strength A for the 
transition 4 <-> 5 (circles) and 5 <-» 6 (squares). Lines are a guide for the eye. (d) 
Same as in (c) but for the transition 6 <-» 7. The insets in panels (c,d) represent 
the structure of the dot GS WF in the strongly correlated regime. Conductance 
unit g = 1.6 ■ 10- 7 r s/GQ 



where the tunnelling electron should enter in order to provide an optimal matching 
between the dot states and obtain a good transmission through the dot. On the other 
hand, the most likely position of the tunnelling electron is essentially dictated by |AL|, 
as the QPWF shows. As a result, a higher conductance is obtained in situations where 
the QPWF is peaked so as to provide a maximal overlap of the dot WFs. With these 
considerations, let us now reexamine figure |5[c,d). 

For the transition 4 <-> 5, as A increases, the dot WFs build up radial and subsequently 
angular correlations, ending up eventually in a molecular state with a square (N = 4) 
or pentagon (JV = 5) symmetry with always a ring-like structure. As such, maximum 
overlap is achieved when the tunnelling electron jumps to the edge of the dot. This is 
the case, in agreement with the results of the QPWF, as confirmed by figure[7Ja). 
The transition 5^6 displays a more interesting double-slope feature. For small A, 
both the dot GS WFs for N = 5 and TV = 6 display weak correlations and have a 
ring-like structure. Similar to the case discussed above the tunnelling electron, enter- 
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Figure 7. Radial behaviour of the modulus of the QPWF |v>(r)| - see II15I I - 
for A = 1.2 (red), A = 1.8 (green) and A = 2.4 (blue) and different transitions 
between dot GSs: (a) 4 <-> 5; (b) 5 <-> 6; (c) 6 <-> 7. At the top of the panels, |AL| 
for the given transition is quoted, see table [T] 



ing at the edge of the dot, provides an optimal overlap of the dot WFs. Therefore, 
a slope similar to the one observed for 4 <-> 5 is obtained for small A. On the other 
hand, for A > 2 one electron is shifted towards the center of the dot. Eventually, the 
WF for N = 6 acquires the shape of a centered pentagon, see the inset in figure [6][c). 
The optimal overlap would be achieved with the tunnelling electron jumping to the 
center of the dot. This however is not allowed for dynamical reasons, as shown by 
the QPWF in figure EJb): the tunnelling electron needs to enter into the dot edge. 
Therefore, an additional suppression of the conductance occurs, which is detected in 
the sharp change of slope of Q shown in figure EIc). 

In the case of 6 <-» 7, for A < 1.8 correlations in both the WFs are weak and the latter 
exhibit a ring- like shape as in all the low A regimes already discussed. For A « 1.8, 
the GS with N = 7 begins to shift one electron towards the center of the dot, while 
the GS of N — 6 remains annular. Since the QPWF is peaked at the edge of the dot, 
this corresponds to a slight, yet noticeable, suppression of Q visible in figure [Hd) for 
A = 1.8. For A > 2, also the GS for N = 6 has one electron in the center of the dot. 
The optimum overlap is again achieved for tunnelling at the dot edge and therefore 
one could expect a return of a power law similar to the one observed for A < 1.6. On 
the contrary, one observes a steeper slope. In order to explain this phenomenon we 
need to consider in detail the spin structure of the dot GSs. 

We consider here the two states \N,L,S,S Z ) = |6, 0,0,0) and |7, 2, 1/2, 1/2). Figured] 
shows the spin-resolved two body correlation function pl z ' s *(r,r') in (|20|) . calculated 
along the outer ring of the quantum dot GSs for TV = 6, 7 and a representative value 
A = 2.4. One electron with spin down is assumed to lie at r' = ro(N, A) and 9' = 0. 
Starting with N = 6, panel (a), one can see that the probability of finding another 
spin-down electron (p^H®)) 1S peaked around 9 w 47r/5 and 9 ps Qtt/5, while the 
one for an electron with spin up (pl ,!t (9)) is highest at 9 = 2ir/5 and 9 = 8ir/5, 
although two relative maxima for this probability are observed also for 9 = Air /b and 
9 = 67r/5. For N = 7 one finds that the probability of finding another spin-down 
electron (P2^(9)) is largest at 9 = 2ir/3 and 9 = 4ir/3, while spin-up electrons are 
maximally likely at 9 = 7r/3, 7r, bir/3. 

In the Wigner molecule regime, only a few different spin configurations may contribute 
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Figure 8. Angular dependence of the spin-resolved two-body density functions 

p s 2 *' s ' z {0) = pl°" s '*(r o (N,\),0,r o (N,\),0) (units: i~ A ) with s' z = -1/2 and 
s z = —1/2 (red line) or s z = +1/2 (green line) for the dot states (a) |6, 0,0,0); 
(b) 1 7, 2, 1/2, 1/2). Here, A = 2.4. 




I t i 




Figure 9. The possible spin pattens, most relevant for linear transport, 
contributing to the GS WF for A = 2.4 and (a,b) AT = 6; (c) N = 7. The 
probe electron lies at 8 = and the white dot denotes the preferred site for 
tunnelling. 



to the dot WF. In the case of N = 6, once the spin direction for the electron at the 
center of the dot is chosen, only two possible spin arrangements are possible for the 
pentagon at the edge. These are shown in Figs. [9^a,b) for the case of a spin up in the 
dot center. Since the WF for N = 6 is a spin singlet, also the two other configura- 
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tions, obtained flipping all spins in panels (a) and (b) are possible (not shown). Since 
the correlation function in figure [S^a) for parallel spin-down electrons is more peaked 
around 9 = 47r/5 and 9 = 6ir/5, one can anticipate that the configuration represented 
in panel (b) contributes more than the one in panel (a). Also for TV = 7 several spin 
configurations for the dot edge exist, once the spin in the center of the dot has been 
fixed. However, the clear peak structure of figure[8jb) strongly suggests a well-defined 
texture of alternating spins in the outer ring of the molecule with a corresponding 
spin-up electron in the center of the dot, consistent with S z = 1/2. Therefore, we can 
infer that the dot WF for N = 7 has the spin structure shown in figure[S{c). Note that 
the discussion for the other states of the multiplet for N = 7 is identical, provided 
that one flips all spins for the states with S z = —1/2. 

Let us go back to the transport properties of the dot. Among all the possible spin 
configurations for jV = 6, the ones with a spin down at the dot center, obtained by 
flipping the spins of those shown in figure [5Ja,b), provide a very poor overlap with 
the state with N — 7 and therefore can be neglected. The configuration shown in 
figure [UJa) also provides a negligible overlap: there is no position around the edge 
for the tunnelling electron so that the final state has the same spin pattern shown in 
figure [9jc) . Concerning the situation shown in Fig [9jb) , the tunnelling electron can 
only jump in the proximity of 9 — tt (white dot in the figure) since all other positions 
would lead to a wrong spin pattern on the edge. This results in a suppression of the 
tunnelling amplitude as compared to the case of small A < 1.8, when the tunnelling 
electron is free to delocalize around the ring due to the negligible correlations of the 
dot WF. 

It is important to note that the situation described before does not occur either for 
4^5, whose conductance is featureless, or for 5 <-> 6 whose change of slope is mainly 
related to the localization of one electron in the center of the dot. 

From the above discussions one can conclude that the transition towards the Wigner 
molecule, accompanied by qualitative rearrangements of the charge or spin textures 
of the dot WF, may be detectable in the linear transport properties. This seems par- 
ticularly relevant when transport involves states with higher numbers of electrons and 
intricate spin patterns, such as N = 6, 7, due to the complex internal structure. Sim- 
pler configurations such as the ones for the transition 4 <-> 5 discussed above may not 
cause any signature in transport. It is worth to notice that the spin effects discussed 
above are subtler than the more common type-II spin blockade [73|,l74j. In the latter, 
the current flow is blocked due to the impossibility to fulfil total spin conservation 
by tunneling events. In the Wigner molecule regime, on the other hand, even if spin 
conservation is satisfied an additional suppression of the current as A increases occurs, 
due to the peculiar internal spin structure of the dot WFs. 



3.2.2. Nonlinear transport In the nonlinear regime, transport also triggers the 
population of excited states of the quantum dot. In this section we discuss one 
particular case, to show signatures of the transition towards the Wigner molecule. To 
be specific, we will concentrate on the regime where only states with A = 5, 6 electrons 
in the dot are involved. Furthermore, we assume strong relaxation: H^>ro, see (|13l) . 

The numerically evaluated Q as a function of V and V g is shown in figure [TUT a) for 
A = 1.4. It exhibits conductance lines corresponding to transitions between the dot 
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Figure 10. (a) Nonlinear conductance Q (units Go same as in figure [6jl as a 
function of applied voltage V (units mV) and gate voltage l/ g (units mV) for 
A = 1.4, T = 750 mK and W = 50Fq. (b) Scheme of the transitions involving 
the ground states for N = 5, 6 and the lowest lying excited states for N = 6. The 
dashed line represents a spin-blockaded transition not visible in the sequential 
regime. 



Table 2. States and corresponding quantum numbers involved in transport 
dynamics for the range of V and V g considered in figure ITOl 



State 


N L 


S 


S z 


\A) 


5 ±1 


1/2 


±1/2 


|B) 


6 








|C> 


6 ±1 


1 


0,±1 


\V) 


6 


2 


0,±1,±2 



GSs or between GS and lowest-lying excited states. The scheme of the expected 
lines is shown in figure fTOT b) for the voltages region considered here. The transitions 
corresponding to each line are shown, with their quantum numbers given in table O 
The blue lines represent transitions between the dot GSs for TV = 5 (|-4)) and TV = 6 
(\B)). The red lines represent channels involving the GS of TV = 5 (|-4)) and one of 
the first two excited multiplets of TV = 6: the lowest one is denoted as \C) while the 
next-to-lowest is \T>). Since calculations are performed for temperatures smaller than 
the average level spacing between the dot multiplets, in the strong relaxation regime 
transitions among the excited states of the dot cannot occur. Each of these transition 
lines corresponds to the opening of the specific transport channel involving an excited 
state of the dot with TV = 6. Note that transitions involving excited states for TV = 5 
are not present in the considered range of V and V g , since they lie at higher energies. 
A comparison between the scheme of figure HUT b) and the calculated conductance, 
figure llOf a) shows that only the first transition line (red solid), corresponding to 
\A) — > \C) is observed, while the one corresponding to the second excited multiplet of 
TV = 6 (red dashed) is absent. By inspecting table one notes that the transition 
\A) — ► \T>) involves |A5| > 1/2 and therefore is forbidden [73l [74], leading to a 
vanishing conductance. 

Figure [TT] shows plots of Q as a function of the applied voltage for different values of A. 
In all panels, the value of Vg has been chosen to lie between the green square and the 
green dot in figure fTOT b). The peak at lower V in Figs. [PIT a-d) corresponds to the GS 
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Figure 11. Nonlinear conductance as a function of the applied voltage V (units 
mV) near the transition 5 <-> 6 for (a) A = 1.4, V g = 39.3 mV, T = 750 mK; 
(b) A = 1.8, V g = 27.3 mV, T = 150 mK; (c) A = 2.2, V g = 20.2 mV, T = 150 
mK; (d) A = 2.6, V g = 15.9 mV, T = 100 mK. In all plots, W = 50T and the 
conductance unit Qq is the same as in figure [6] 



to GS transition |_4) — > \B), while the second one to the transition \A) — ► \C). As is 
visible from the voltage ranges of the plots, the dot level spacing (corresponding to the 
distance between the nonlinear conductance peaks) gets narrower as A is increased. 
In order to be able to resolve both conductance peaks, calculations for higher A have 
been performed at lower temperatures than those at smaller A. 

Comparing the panels (a) and (b) for A < 2 with panels at A > 2 (c) and (d), 
a qualitative difference in the behaviour is easily observed: for weaker interaction 
strengths, the first peak is always higher than the second one, while for A > 2 the 
situation reverses drastically. 

Such a behaviour cannot be attributed to the difference in temperature between 
different calculations. Indeed, calculations for A = 1.4 and A = 1.8 performed at 
lower temperatures display narrower conductance peaks but still with almost equal 
height. Increasing T for A = 2.2 and A = 2.6 always suppresses the conductance for 
the GS to GS peak with respect to the transition towards the excited state in all the 
temperature range in which the two peaks are resolved. 

The behaviour of the nonlinear conductance can be related to qualitative changes 
in the WFs of the excited states for N = 6. Figure [12] shows the states involved in 
the transport dynamics of the two conductance peaks discussed above. The height of 
each peak is determined by the available transport channels and by the transition rates 
connecting the dot states. As the interaction strength is increased the first conductance 
peak, involving only transitions between GSs - see figure [T27 a) - behaves exactly as 
the linear conductance peak discussed in Sec. 13.2.11 The height of the second peak is 
on the other hand determined by two families of transport channels connecting the GS 
with TV = 5 to the excited multiplet of N = 6: as shown in figure fTSTb). channels with 
either AL = or |AL| = 2 are possible. The corresponding transition amplitudes are 
shown in figure [T3T a) as a function of A. Transition rates are proportional to these 
amplitudes, see (fl"0|) . One can see that for A < 2 the transition channel with AL = 
is strongly suppressed, while the one with |AL| = 2 is larger and decaying with A. For 
A > 2, a sudden decrease of the transition amplitude for the channel with |AL| = 2 is 
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Figure 12. (a) Scheme of the transitions involved in the first conductance peak 
shown in figure [TT] (b) Same as in (a) but for the second conductance peak. 



found, while the one for the channel with AL = jumps to a very large value. 
This peculiar behaviour can be again explained using the QPWFs. Figure [TH shows 



\Oi\ 



□ AL = 
olALI = 2 
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Figure 13. (a) Squared modulus of the transition amplitude |d| 2 - see iTTTJ - 
for the processes |5, ±1, 1/2, ±1/2) -> |6, =pl, 1, ±1) with |AL| = 2 (circles) and 
|5, ±1,1/2, ±1/2) -* |6,±1,1,±1) with AL = (squares) shown in figure [T2TbV 
The insets show the character of the dot WF for high values of A. (b) Same 
as in (a) but for the GS to GS transition |5, ±1, 1/2, ±1/2) -> |6, 0,0,0) shown 
in figure 1X27 a). Amplitudes for transitions involving triplet states with S z = 
are proportional to the ones shown here, the proportionality constant being a 
Clebsch-Gordan factor l/\/2. Amplitudes for transitions with |AS^| > 1/2 are 
zero. 



the QPWF for the transition with AL = (panel a) and for the one with AL = 2 
(panel b). When A < 2, the case of AL = has nonzero amplitude near the center of 
the dot and consequently a very small overlap between the configurations for TV = 5 
and N — 6 which have a ring-like character. This results in a very small transition 
amplitude. On the other hand, the case of |AL| = 2 has a large QPWF near the edge 
of the dot and therefore provides a much better overlap between the dot configurations. 
The situation reverses dramatically for A > 2: the developing radial correlations induce 
a qualitative change in the dot WF for TV = 6 and one electron moves near the center 
of the dot, in analogy to the case of the GS. In this case, a much better overlap occurs 
for AL = 0, since the extra electron preferably sits in the center of the dot, see the 
right insets in figure HUT a). The transition with |AL| = 2 is strongly suppressed since 
the corresponding QPWF is cx r 2 and thus negligible near the dot center. 
Therefore, for A < 2 transport through the excited state occurs essentially via the 
channel with |AL| = 2, whose amplitude is similar to the one for the GS to GS 
transition, see figure [T3jb). This explains why, for A < 2, the two conductance peaks 
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Figure 14. Radial behaviour of the modulus of the QPWF \<p(r)\ for A = 1.2 
(red) A = 1.8 (green) and A = 2.4 (blue) and different processes: (a) 5 <-> 6 GS to 
first excited with AL = 0, i.e. |5, ±1, 1/2, ±1/2) <-» |6, ±1, 1, ±1); (b) 5 <-» 6 GS 
to first excited with | AL\ = 2, i.e. |5, ±1, 1/2, ±1/2) <-> |6, =f 1, 1, ±1). Note that 
in (a), the red and green curves for A = 1.2 and A = 1.8 have been magnified by 
a factor 20. 



in figure ITTT a.b) have almost equal height. For A > 2 the channel with AL = clearly 
dominates, due the peculiar rearrangement of the dot WF for N = 6 moving towards 
the Wigner molecule. The amplitude for this channel is larger than the one for the 
GS to GS transition and therefore the conductance peak for the transition involving 
the excited state is higher than the one for the GS to GS transition as shown in 
figure [TTTc.d). 

4. Conclusions 

In this paper we have investigated correlation effects in a quantum dot via linear 
and nonlinear transport, employing the PHF technique. For increasing interaction 
strength, the ground and excited dot states have been analyzed for 4 < N < 7. 
As the strength of the Coulomb interactions increases, the dot WFs build up radial 
and angular correlations, smoothly crossing over from a liquid-like regime to Wigner 
molecular states. Most strikingly, we have demonstrated that signatures of such a 
crossover may appear both in the linear and in the nonlinear transport properties. 
These signatures have been interpreted with the systematic study of both two-body 
correlation functions and QPWFs. 

In the linear regime, we have observed an exponential suppression of the conductance 
as the transition towards the Wigner molecule takes place. In cases when the latter is 
accompanied by strong qualitative rearrangements of the dot WFs, strong mismatches 
of the dot WFs involved in the transport process may occur. This leads to a stronger 
suppression of the conductance, as observed for the case 5 <-> 6. A mismatch of the 
dot WFs due to the emergence of particular spin structure of the dot states may 
also occur, as exemplified by the case of 6 <-*■ 7. Also this fact leads to an increased 
suppression of the linear conductance. 

In the nonlinear regime, the conductance may even be enhanced by the formation of 
a Wigner molecule within the dot, as shown by the study of the transport dynamics 
of the lowest-lying excited states for N = 6. 

The effects described above are due to the qualitative rearrangements of the charge 
or spin patterns of the dot states, occurring during the transition towards the Wigner 
molecule regime. As a possible extension of this investigation, it would be interesting 
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to devise a method to investigate the internal spin structure of the Wigner molecule, 
e.g. by analyzing the effects of spin-dependent tunnel barriers. Effects due to Wigner 
molecules should have a profound impact in coherent regimes and could lead to strong 
signatures detected by analyzing e.g. the cotunnclling regime. Finally, it would be 
interesting to consider the effects of applied magnetic fields, which are known to 
strongly modify the properties of Wigner molecules. 

We expect that results similar to the ones shown in this paper hold also for planar 
quantum dots and that they could be in principle observed experimentally. 
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Appendix A. Tunnelling Hamiltonian 



In this appendix we derive H t for a quantum dot connected to two external emitter or 
collector leads via tunnelling barriers located around z = z a (a = E, C), see figure [1] 
Our starting point is the Bardeen tunnelling Hamiltonian |83] 



i 

2m* 



a=E,C 



dRdz S( z — z n ) 



OZ 



h.c. 



(A.l) 



Here, R = (x, y) and z a is a point within the tunnelling barrier between the lead a 
and the dot. Furthermore, ^(R, z) is the system field operator expanded into emitter, 
dot and collector contributions as ^(R, z) = >]/d(R, z) + J2 a =E c ^aCRj z ) with 

* Q (R,z) = 5>£>(R,z) 
*d(R, z) 



E 

V 



^\R,z)d v . 



(A.2) 



(A.3) 



Here, $^ Q ^ (R, z) and (R, z) are a set of single particle eigenfunctions for leads and 
dot, respectively. Note that these sets may be complete but leads and dot states need 
not be orthogonal. In the following, we assume that the longitudinal (z) and transverse 
(x, y) motions are decoupled, which allows the factorizations of the single particle 
WFs. Furthermore, we will concentrate on a system with cylindrical symmetry. For 
the leads one has $^(R, z) = (f> Ua (R)xfc„ ( z ) where k a is the momentum along the z 
direction and v a is a set of quantum numbers describing the transverse motion. For 



the quantum dot we choose ^^(R, z) = ^(R)xd(z), where 77 



.CD) 
5zr} 



} 



collectively denotes the FD quantum numbers. Choosing the basis of FD states is 
not restrictive: as shown in |Appendix C[ every orthonormal and complete basis for 
the single particle states of the dot is equivalent. We assume a sufficiently tight 
confinement in the z direction so that longitudinal motion of electrons in the dot is 
effectively frozen into the lowest subband with WF xd(-z). The longitudinal dot WF 
Xd{z) is evanescent within both tunnelling barriers. For the leads, Xk E (z) (Xk c ( z )) 
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is evanescent under the tunnelling barrier near the emitter (collector) and essentially 
zero near the collector (emitter). Substituting the explicit expressions of the field 
operators into (|A.1[) and taking into account the form of the single particle WFs one 



obtains H t = H, 



(0) 



h[ 1] with 



H 



(0) 



H 



(i) 



= E E*&AeAft 

a=E,C 

= E Efi4,e^ 



h.c. 



(A.4) 
(A.5) 



where 



A (D) 



J") 



2m" 
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2m* 



2m* 



E 

a=E,C 



dz 

dxD(z) 



dz 



- Xk' a {z a ) 
Xd(zq) 



dxl> (z) 



dz 

dxh(z) 



dz 



Xl a (z a ) 



dxv(z) 



dz 



Xo(z a ) 



dxt(z) 



dz 



f dR ^JR)^(R) 



(A.6) 



The first term in H^ produces a weak, one-body scattering within the leads, while 
the second term gives rise to a small uniform shift of the dot energy levels. Both 
effects can be safely neglected. The relevant term is H\ f , which produces scattering 
of electrons between the leads and the dot. For a pillar quantum dot, assuming 
harmonic confinement of electrons in the emitter and collector, with frequency uj a , 
the WFs </v a (R) are FD states. We consider the case ui a « w. Equation (|A.6[) yields 
essentially 



- 4 a h 



(A.7) 

where in the last equality we have neglected the weak dependence of the tunneling 
matrix element on k a . Finally, the tunnelling Hamiltonian assumes the form 



(A., 



Appendix B. Tunnelling rates 

In the sequential regime, tunnelling rates between initial \T) and final \T) states of 
the system are obtained via the Fermi's golden rule [84] 

Tx^r = 2tt |(^|i? t |X)| 2 5{E? - Ex) , (B.l) 
where Ex {Ejr) is the total energy of the system in the initial (final) state and 
\X) = |2e) <8> \Xc) <8> |2b) ; \F) = <8> \Fc) <8 \Tu) , 

where \I a ) ( \T a )) is the initial (final) state for lead a while \It>) (I-^p)) is the dot 
initial (final) state, with electron number Nx D (A^ D ) and energy £?x D {E^ D ). Note 
that dot states can be either ground or excited. In order to have a nonvanishing 



Transport properties of quantum dots in the Wigner molecule regime 



24 



contribution in the sequential regime, \Njr D — Nj D \ = 1 must hold. This implies 
that (|B.1[) is diagonal in the barrier index a and that sequential tunnelling events 
through the barriers are independent. The tunnelling rate has the general structure 

rw= E E r ^, 

a=E,Cp=±l 

where the contribution with p = +1 (p = — 1) represents tunnelling into (out from) 
the dot via lead a. 

Since we are interested into the dot dynamics only, we perform a thermal average over 
\T a ) and a summation over \T a ) obtaining transition rates among dot states only 

Txo^o = E E rt^o ■ (B-2) 

a=E,Cp=±l 

The leads are assumed to be in equilibrium with respect to their electrochemical 
potentials fj, a = fio + <5/x Q , where <5/x Q is a shift due to the presence of an applied bias 
voltage V. In the case of a pillar dot, see (|A.8j) . one obtains 



2d — >^ 7 D 



r (Q) \o p \ 2 / p (/i D - fi a ) , 



(B.3) 



where O v 



and / p (.E) = p/(S) + (l-p)/2, with f(E) = [l + cxp^S)]" 1 the Fermi distribution at 
inverse temperature (3 = (fee the Boltzmann constant) and /id = Ejr — Ej D . 

Appendix C. Equivalence of different bases 



In this appendix we will show the equivalence of different single particle bases of the 
dot for the definition of the tunneling Hamiltonian. Consider two orthonormal and 
complete bases of single particle dot states labeled by {rf\ and {/x}, with corresponding 
orbitals 0^(R), (/>^(R) and Fermi operators d v , (in this section we omit over hats 
from the operators to simplify the notation). They are connected by a unitary matrix 
M, 



'ii' 



,(R) = E M w^(R) ; ^ = E m w j m- 



(C.l) 



The tunneling Hamiltonian, expressed in the basis {/x} reads 



(i) 



with 



=(«) 



E Et } A,A 



a=E,C f Q ,fi 



2m* 



XL (Za) 



dz 



Xd(^q) 



dxl a {z) 



dz 



dR # q (R)^(R), 
see (|A.6I) . Inverting (|C.1[) for d^ and plugging into (|C.2[) one obtains 



Hi 



(i) 



E E 



5> 



,==(«) 



h.c. 



(C.2) 



(C.3) 



(C4) 
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By virtue of (|C.ip . one can easily see that 

which shows that ()C.2[) is identical to (|A.5j) . This implies that all results are 
independent of the choice of the single particle states for the dot. 
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